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Abstract. We present spectroscopic observations obtained with the FORS1 instrument on the ESO VLT under good seeing 
conditions of the radio-quiet 2 em = 3.04 quasar Q1205-30 and its associated extended Lya emission. The extended Lya 
emission was originally found in a deep narrow band image targeting a z a bs ~ z em Lyman-limit system in the spectrum of the 
QSO. Using spectral point-spread function fitting to subtract the QSO spectrum, we clearly detect the extended Lya emission 
as well as two foreground galaxies at small impact parameters (2.12 ± 0.04 and 2.77 ± 0.07 arcsec). The redshifts of the two 
foreground galaxies are found to be z = 0.4732 and z = 0.865. We determine the redshift and velocity profile for the extended 
Lya emission, and analyzing the velocity offsets between eight QSO emission lines we refine the quasar redshift determination. 

We use the new redshifts to infer the geometry of the complex. We find that the extended Lya emission is clearly associated 
with the quasar. A Lya luminosity of 5.6 x 10 43 erg s _1 places this extended emission at the high luminosity end of the few 
previous detections around radio-quiet quasars. The extended Lya emission is best explained by hydrogen falling into the dark 
matter halo inhabited by the quasar. 
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1. Introduction 

Quasar host galaxies are visible tracers of the close environ¬ 
ment of this powerful type of active galactic nuclei (AGN). 
The feeding of the central engine from the host galaxy, and the 
feedback of the quasar to the host are important unknown fac¬ 
tors in current numerical models, which need to be understood. 
The study of quasar host galaxies is difficult, because of the 
high contrast between the bright point-source quasar and the 
faint, extended host galaxy. Surveys of AGN host galaxies have 
primarily targeted radio-loud quasars (RLQs) and radio galax¬ 
ies (RGs) (e.g. Lehnert et al. 1992, 1999; Reuland et al. 2003; 
Sanchez & Gonzalez-Serrano 2003), despite the fact that the 
majority of quasars are radio-weak or radio-quiet. In such sur¬ 
veys it has been found that RLQs and RGs at low z (here taken 
to mean z < 1) reside in luminous elliptical galaxies, while 
radio-quiet quasars (RQQs) are found in both elliptical and 
early spiral galaxies. At intermediate redshifts (1 < z < 2) 
galaxies in general appear to have more disturbed morpholo¬ 
gies, making it difficult to apply the simple classification of 
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“ellipticals” and “spirals”. Studies have found that host galax¬ 
ies of RQQs are 1 — 2 mag fainter than hosts of RLQs with 
similar luminosity (Falomo et al. 2001; Kukula et al. 2001). 
Similarly, in a study of quasars out to 2 = 2.1 the black holes 
of RLQs are typically found to be 45% more massive than their 
radio-quiet counterparts (McLure & Jarvis 2004). At 2 > 2 
the cosmological surface brightness dimming makes it increas¬ 
ingly difficult to make secure detections of the host galaxies. 
For RLQs some 20 examples of 2 > 2 hosts are seen (Hu et 
al. 1991; Heckman et al. 1991a,b; Steidel et al. 1991; Lehnert et 
al. 1992, 1999; Wilman et al. 2000), but only few surveys have 
targeted RQQs and with limited success (Bremer et al. 1992; 
Lowenthal et al. 1995; Fynbo et al. 2000a; Mpller et al. 2000; 
Ridgway et al. 2001; Bunker et al. 2003). The relatively faint 
hosts of RQQs compared to RLQs makes them more difficult 
to detect. A promising method is narrow-band imaging tuned 
to the Lya line at the quasar redshift (Hu & Cowie 1987; Hu 
et al. 1996). Haiman & Rees (2001) predict that gas enshroud¬ 
ing a quasar between redshifts 3 and 8 would be photoionized 
by the quasar UV emission and should be detectable at a sur¬ 
face brightness of 10 -18 to 10 -1 ' erg s _1 cm -2 arcsec -2 in 
the Lya line. These limits have only been reached for very few 
surveys. 
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Table 1 . Observation log. 



PA 

Grism 

Seeing 

Exp-time 

PA1 

67.20° 

G600B 

0'.'91 

4 x 1800 s 

PA1 

67.20° 

G600R 

0769 

2 x 1800 s 

PA2 

7.90° 

G600B 

0769 

4 x 1800 s 

PA2 

7.90° 

G600R 

07 73 

3 x 1800 s 


In this paper we report on a spectroscopic study of the 
sightline towards the radio-quiet quasar Q1205-30 at z = 3.04 
and its associated extended Lya emission detected by Fynbo et 
al. (2000b) (hereafter Paper I). 

The paper is organized as follows: In Sect.|2]we present the 
observations and data reductions. We continue in Sect. 0with 
the results of our analysis. In Sect. ^ we discuss the foreground 
galaxies, and we present a simple model for the extended Lya: 
emission in Sect.|3 We end in Sect. [6] with a discussion of the 
origin of the extended emission. Unless stated otherwise we 
will use n m = 0.3, T2 a = 0.7, 7/o = lOO/i km s _1 Mpc _1 = 
70 km s _1 Mpc _1 . In this model a redshift of 3.04 corresponds 
to a luminosity distance Dl = 25.8 Gpc and a distance mod¬ 
ulus of 47.1. One arcsec on the sky corresponds to a projected 
distance of 7.67 proper kpc and the lookback time is 11.4 Gyr 
(84.5% of the time since Big Bang). 


2. Observations and data reduction 

The observations were carried out with the Unit Telescope 1 
(Antu) of the ESO Very Large Telescope (VLT) on March 
4-5, 2000, under photometric and good seeing conditions. 
The data were acquired with the FOcal Reducer/low dis¬ 
persion Spectrograph (FORS1) instrument in Multi Object 
Spectroscopy (MOS) mode with the red G600R and blue 
G600B grisms as a part of a larger campaign (see Fynbo et 
al. 2001). The slitlet used was 1.2 arcsec wide and 20 arc- 
sec long, and the positions on the CCD resulted in a wave¬ 
length coverage of approximately 4000 — 6000 A (G600B) 
and 5800 — 7800 A (G600R). The slitlets had position an¬ 
gles of 67.20° (PA1) and 7.90° (PA2) East of North centred 
on the quasar (see below). For the observations we used the 
standard resolution collimator. During observations with the 
G600B grism the CCD was binned 2x2, and resulting in a 
pixel size of 0.4 arcsec by 2.4 A for G600B and 0.2 arcsec by 
1.2 A for G600R. The seeing in the combined science frames 
was measured to be 0.9 arcsec (PA1) and 0.7 arcsec (PA2) at a 
wavelength of ~ 5100 A, leading to spectral resolutions of 5 A 
(PA1) and 4 A (PA2). The exposure times are given in Tabled 
Contour plots of the 18 x 18 arcsec 2 field of Q1205-30 
imaged in I and Lya narrow band is shown in Fig. ^ taken 
from Paper I. We shall here follow the naming convention of 
that paper, i.e. gl is the blue galaxy SW of the quasar, g2 is 
the red galaxy NE of the quasar, and S6 is the extended Lya 
emission N and NE of the quasar. The projected distances from 
the QSO are 2.12 ± 0.04 arcsec for gl and 2.77 ± 0.07 arcsec 
for g2. For PA1 the slitlet covers gl, g2, the QSO and a part of 
S6, whereas for PA2 the slitlet covers the central part of S6 and 
the QSO. 





Fig. 1 . The 18 x 18 arcsec 2 field of Q1205-30. Upper panels: 
Contour plot of the /-band image before (a) and after (b) PSF 
subtraction of the quasar. North is up and East is left. Lower 
panels: Contour plot of the narrow band image before (a) and 
after (b) PSF subtraction of the quasar. North is up and East is 
left. Taken from Paper I. 


2.1. Basic reductions 

The individual science frames were bias subtracted using stan¬ 
dard techniques. The flat fielding was done by first filtering the 
flat fields along the dispersion axis with a 61 x 1 pixels (146 A 
long) median filter for the G600B grism, and a 121 x 1 pixels 
(145 A long) median filter for the G600R grism. Then the flat 
fields were normalized by dividing the unfiltered flat fields by 
the filtered ones, and finally we divided the science frames by 
these normalized flat fields. In order to obtain a mean sky spec¬ 
trum spatial bins on both sides of the quasar spectrum were 
filtered using a 1 x 13 pixels median filter (i.e. only filter¬ 
ing along the spatial axis) to remove cosmic ray hits and av¬ 
eraged. Regions used for determining the sky spectrum were 
never closer than 5 arcsec to the QSO on the side of the ex¬ 
tended emission and 4 arcsec on the opposite side. The mean 
sky spectrum was expanded to a two dimensional spectrum by 
duplicating the ID spectrum and subtracted from the unfiltered 
science frame. 

2.2. Spectral extractions 

The science frames were coadded and the quasar spectrum was 
optimally extracted using the code described in Mpller (2000). 
After spectral point-spread-function (SPSF) fitting and removal 
of the QSO, the extended Lya emission of S6 was clearly visi¬ 
ble at both PA1 and PA2 (see Fig. [3. 

PA1 was aligned with the two galaxies gl and g2 (see 
Fig. [D and centred on the QSO, so in this case we had to de- 
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Fig. 2. Spectra at PA1 (left column) and PA2 (right column), a) - b) The 2D quasar spectrum, c) - d) The 2D quasar spectrum 
has been subtracted between 4850 A and 5000 A revealing the extended Lya emission. The zero point of the spatial axis is set 
to the centroid of the QSO. Note the large residuals from the quasar SPSF subtraction, which are due to the large shot noise near 
the centroid, e) - f) The spectrum of the QSO Lya emission line (light grey curve) and the spatially averaged spectrum of the 
extended Lya emission (black curve). The maximum flux of the extended Lya emission is scaled to 90% of the maximum QSO 
flux. The total line flux in extended Lya is around 1% of the flux from the QSO Lya line. The overall flux scale in this plot is 
arbitrary. 


compose the spectrum into its individual components. We em¬ 
ployed an iterative procedure to separate the contributions from 
gl and S6 to the total quasar flux: i) Extract and remove the 
QSO spectrum, ii) Extract and remove the spectrum of gl and 
S6. The procedure converged to a stable solution after three it¬ 
erations. For the PA1 observations using grism G600R the see¬ 
ing conditions were slightly better than for those using G600B 
(see Table Q, so the projected distances of gl and g2 were 
enough to bring them outside the QSO point-spread function 
(PSF) and we could make a normal extraction of the quasar 
spectrum. 

In the PA2 spectra we had only the QSO and S6 on the 
slit. In order to make sure that the extended flux of S6 was 
not modifying the quasar SPSF, we here used an option in the 
code which allowed us to exclude a wavelength region from the 
construction of the quasar SPSF. In the region from 4902 A to 
4987 A the code therefore only fitted and extracted the QSO 
spectrum, it did not update the SPSF lookup table (for details 
see Mpller 2000). For the general problem of decomposing a 
2D spectrum of two superimposed objects there is a certain de¬ 
generacy of solutions. One may choose to assign the maximum 
amount of flux to one object, to the other object, or to aim for 
somewhere in between. In this case we decided to opt for a so¬ 


lution in between, which avoids digging a hole in the extended 
Lya emission near the QSO centroid. The degeneracy has no 
consequence at distances larger than 1 arcsec from the QSO, 
but closer to the QSO the Lya surface brightness of S6 is very 
uncertain. 

The output from the code is the optimally extracted ID and 
2D spectrum of the QSO as well as the 2D spectrum of the 
extended Lya emission. The 2D spectrum of the extended Lya 
emission and the ID spectrum of the QSO (with cosmic ray hits 
removed) are shown in Fig. Eh® and Fig.0 respectively. For 
comparison we show a zoom of the Lya emission of the quasar 
and the spatially averaged extended Lya emission in Fig. ctcr. 

2.3. Wavelength and flux calibrations 

The spectra in Fig. 0 we re wavelength calibrated using the dis- 
pcor task in IRAF 1 . The RMS of the deviations from a 4. th 
order Chebychev polynomial fit to 12-19 lines were 0.6 A for 


1 IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of Universities 
for Research in Astronomy, Inc., under cooperative agreement with 
the National Science Foundation. 
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Fig. 3. Extracted spectrum of Q1205-30. The most prominent emission lines are indicated. The spectrum has not been extinction 
corrected nor flux calibrated. 
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Fig. 4. Surface brightness and velocity profiles of the extended Lyct emission based on spectra at PA1 (left column) and PA2 (right 
column), a) - b) Surface brightness profile. The error bars are based on photon statistics only. Systematic errors due to the flux 
calibration are not considered. The thick dashed line shows our 5cr detection limit, c) - d) Velocity profile of the extended Lyct 
emission measured relative to the systemic redshift of the quasar (see Sect. 13.2> . The uncertainties are statistical uncertainties on 
the position of the maximum of a fitted Gauss curve (Landman et al. 1982). In both PAs the most distant detection was too faint 
to allow a secure determination of the velocity. A rest-frame velocity of 600 km s~ 1 relative to the systemic quasar redshift (see 
Sect. 13.2> corresponds to a redshift z = 3.049. 


the G600B spectra and 0.08 — 0.16 A for 28-33 lines in the 
G600R spectra. 

The flux calibration was done as follows. First we estimated 
the continuum of the QSO away from the emission lines. We 
divided the spectrum by the continuum, obtaining a flat spec¬ 
trum, and forced it onto an arbitrary power-law. Measuring the 


Bessel B — I colour on the resulting spectrum and compar¬ 
ing to the observed B — I colour (Paper I) we calculated the 
slope of the power-law which would bring the two in agree¬ 
ment (a = 0.65). We then forced the spectrum onto this power 
and normalized it to agree with our photometric measurements. 
The advantage of this procedure is that to the first several or- 
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Table 2. Emission lines of gl. 


Line 

Avac 

Zem 


(A) 


Oil A3727 

5492.21 

0.4733(5) 

H/3 A4862 

7163.67 

0.4731(9) 

Oiii A4960 

7308.07 

0.4733(2) 

Oiii A5007 

7377.66 

0.4731(1) 


ders all absolute and differential slit losses as well as atmo¬ 
spheric absorption are automatically taken into account for all 
point source objects on the slit. We estimate the absolute flux 
calibration to be correct to within 10% while the relative is bet¬ 
ter. This applies to point sources, but not to extended sources 
for which there are additional slit losses. We present surface 
brightness profiles and rest-frame velocity curves of the ex¬ 
tended Lya emission (relative to the systemic redshift of the 
quasar, see Sect. 13. 2> in Fig. 0 

3. Results 

3.1. Redshifts of the galaxies gl and g2 

The primary purpose of this study is to clarify the nature of the 
extended Lya emission at z « 3.04. As pointed out in Paper I 
the galaxies gl and in particular g2 may have a lensing effect, 
enhancing and stretching the Lya patch. Therefore we first set 
out to determine their redshifts. 

We identify four emission lines in the spectrum of the blue 
galaxy gl. The AB-magnitudes of gl are B(AB) = 24.8, 
I(AB) = 23.3, n(AB) = 24.1 (Paper I). Disregarding the 
weak and noisy O III A4960 line, we derive a mean redshift of 
z = 0.4732 ± 0.00011 (Table[2}. In the 2D spectrum (Fig. |5j 
we see along the slit a tilt of the emission lines caused by the 
rotation of gl. The A4960 line is too weak, but for the three 
remaining lines OII A3727, H (3, and O III A5007, we mapped 
out their rotation profiles. The three profiles are identical within 
the errors, so we combined the profiles and performed a joint 
fit (see Fig. [6}. From linear extrapolation of the rotation pro¬ 
file down to the position of the QSO (b = 0 kpc), we find that 
gl would cross the QSO spectrum at —212 ± 20 km s _1 corre¬ 
sponding to z = 0.4722±0.0001. Assuming instead a flattened 
rotation curve after the observed point closest to the QSO, we 
find that gl would cross the QSO spectrum at —100 ± 10 km 
s -1 corresponding to z = 0.4727 ± 0.0001. We shall return to 
a search for absorption at this redshift in Sect. 0] 

The red galaxy g2 has no emission but absorption lines 
(see Fig. |7}. The AB-magnitudes of g2 are B(AB) > 26.5, 
I(AB) = 22.22, n(AB) > 24.8 (Paper I). It was suggested 
in Paper I that g2 most likely is a normal elliptical galaxy at 
^ > 0.5. Minimum-x 2 fitting to redshifted template elliptical 
galaxy spectra (Kinney et al. 1996) gave a best fit redshift of 
z = 0.865 ± 0.003, confirming the prediction of Paper I. The 
observed spectrum and the best-fit template are shown in Fig.Q 
The observed I-band magnitude of g2 corresponds roughly to 
an absolute B-band magnitude of Mb = —22.1, which is 0.7 
mag brighter than M for field galaxies in the redshift interval 
0.75 < ^ < 1.0 (Cross et al. 2004) 




7100 7150 7200 7250 7300 7350 7400 

Observed wavelength (A) 


Fig. 5. a) Spectrum of the blue galaxy gl smoothed with a 
7 A wide boxcar filter. The identified lines are indicated, and 
strong residuals from night sky lines are marked “n-s”. b) 2D 
spectrum of the QSO residuals (at 0 arcsec), the blue galaxy gl 
(at 2.12 arcsec), and the red galaxy g2 (at —2.77 arcsec). The 
rotation profile of gl is clearly visible in the emission lines H (3 
at 7162 A and O III at 7376 A. 


kpc 

6 8 10 12 14 16 



arcsec 


Fig. 6. Rotation profile of galaxy gl mapped out from the 
On A3727 (□), H (3 A4862 (o) and Oiii A5007 (x) lines. The 
straight line is the best combined fit to all velocity profiles. The 
axaxis denotes the distance from the QSO in arcsec (bottom) 
and kpc (top), whereas the //-axis is the velocity relative to the 
redshift of O ill A5007 {z = 0.47311). 


3.2. The systemic redshift of the QSO 

The quasar spectrum contains a Lyman-limit system (LLS) 
very close to the redshift of the quasar (Lanzetta et al. 1991). In 
order to establish the relation between the QSO, the LLS and 
the extended Lya emission we need to determine the precise 
redshifts of all three. For QSOs this is not entirely straight for¬ 
ward. It is well known that high-ionization lines are blueshifted 
with respect to the QSO systemic redshift by typically several 
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hundred km s 1 (e.g. Tytler & Fan 1992). For low-ionization 
lines the blueshifts are known to be small or zero. 

We measure equivalent widths, line fluxes and vacuum- 
corrected wavelengths for eight both high- and low-ionization 
emission lines in the quasar spectrum after removing absorp¬ 
tion lines from intervening systems. The results are presented 
in Tabled Because the LLS redshift is very close to the QSO 
redshift, the quasar Lyo emission line is heavily absorbed, so 
the observed centre wavelength is dependent on how the line 
is reconstructed, resulting in an uncertainty of 20 A. We cal¬ 
culate the systemic redshift of the quasar using the subset of 
the lines listed in Table 0 which also appear with blueshift- 
adjusted rest-frame wavelengths in Tytler & Fan (1992). The 
five lines used are N V, Si II / O I, Si IV / Oiv], C IV, and C III]. 
The inverse-variance weighted average of these lines provides 
a final zqso = 3.041 ± 0.001. 

3.3. Absorption line systems 

In the quasar spectrum we identify 12 C IV absorption systems 
between the Lyo and C IV emission lines. In Fig. [8] we plot 
this section of the QSO spectrum normalized to the continuum. 
The line identifications are presented in Table 0] and the red- 
shifts and identified lines for the 12 systems are summarized 
in Table 0 We expect to detect 5.2 ± 1 ClV absorbers per 
unit redshift (Fig. 5 in Sargent et al. 1988). We find 13.9 sys¬ 
tems per unit redshift between the Lyo (z = 2.17) and Civ 
(z = 3.04) lines of the QSO. Applying the same selection cri¬ 
teria as the S2 sample of Sargent et al., namely rest equivalent 
widths Wrest > 0.15 A for both lines in the CIV doublet, a rest- 
frame velocity relative to the QSO v < —5000 km s 1 (in this 
case corresponding to z a b s < 2.974), and grouping systems 
separated by less than 1000 km s -1 , we are left with 5.8 ± 2.4 
systems per unit redshift. Using these criteria Sargent et al. find 
2.2lg'g systems per unit redshift at 2 = 2 — 3 (Fig. 6 in Sargent 



4000 5000 6000 7000 

Observed wavelength (A) 


Fig. 7. Spectrum of g2 (thin line) smoothed with a 35 A wide 
boxcar filter, and template of an elliptical galaxy (thick line) at 
2 = 0.865 (Kinney et al. 1996). Strong residuals from night 
sky lines are marked “n-s”. 


et al. 1988). Our result is higher than normally observed, but 
still marginally consistent with the Sargent et al. study. 

The LLS at z ~ 3.034, which was the original target 
(Paper I), is identical to the absorption system K, for which 
both low- and high-ionization lines are detected. This is in 
agreement with the findings for other z a b s ~ z em absorbers 
(Savaglio et al. 1994; Hamann 1997; Mpller et al. 1998). By 
fitting line profiles to the Lyo and Lyfj lines we find the HI 
column density to be in the interval 17 < log N\\ , < 19.9. 
The best fit gives log Am = 19.5. In Fig. [3 we plot the 
absorption lines originating from system K in velocity space 
relative to the redshift obtained from the OI absorption line, 
zoi = 3.0331 ± 0.0002. We notice that high-ionization lines 
have systematically higher redshifts than low-ionization lines. 
We use the redshift from the OI line as the redshift of the low- 
ionization region, and we take the redshift obtained from the 
CIV doublet, zciv = 3.0340 ± 0.0002, to be the redshift of the 
high-ionization region. We find that the region of highly ion¬ 
ized elements moves with a velocity of ~ 60 km s -1 relative 
to the low-ionization region. Furthermore, we notice an H I ab¬ 
sorption system in the red wing of both the Lya and Ly [3 lines 
of the LLS (see Fig.|9j. This system, which we will call system 
Kl, has a redshift of 3.039 ± 0.002, and we find no associated 
metal-lines. 

We detect high-ionization N V absorption for the systems 
I and J. If the high degree of ionization is caused by the high 
UV flux from the quasar, the systems must be located between 
the QSO and the LLS at z « 3.034, since no UV photons pass 
through the LLS. The redshifts of systems I and J are lower 
than that of the LLS (system K), which suggests that systems 
I and J are high-velocity clouds. Furthermore, the high degree 
of ionization suggests that they are associated with the QSO. 
The apparent line-locking between N V A1242 of system I and 
N V A1238 of system J strengthens this hypothesis. The line¬ 
locking effect is thought to occur in clouds driven by radiation 
pressure and accelerated via absorption until the wavelength of 
the feature falls in the shadow of another line in a neighbouring 
cloud (Vilkoviskij et al. 1999; Srianand et al. 2002). The veloc¬ 
ities of the absorption systems are —4050 km s -1 (system I) 
and —2850 km s -1 (system J) relative to the systemic redshift 
of the quasar. 

4. Foreground galaxies 

A classical way to search for galaxies at intermediate or high 
redshifts has been to look for absorption systems in the spec¬ 
tra of background QSOs (e.g. Weymann et al. 1979). Studies 
of the galaxy counterparts of Mg II absorption systems have 
concluded that the galaxies responsible for the absorption are 
bright galaxies which can be found at impact parameters of 
~ 100 kpc (Bergeron & Boisse 1991, Guillemin & Bergeron 
1997). Using a sample of damped Lya absorbers (DLAs) 
Lanzetta et al. (1995) found that most galaxy counterparts of 
DLAs at z < 1 have gaseous haloes extending over ~ 230 kpc. 
This is in contrast to studies at high redshift (2 = 2 — 3.5), 
where Mpller et al. (2002) have found that galaxy counter¬ 
parts of DLAs reside at impact parameters ~ 10 kpc. The 
large impact parameters have been used to advocate a picture 
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Table 3. Measured characteristics of eight QSO emission lines. The columns show the line identification, the observed wave¬ 
length in vacuum, the equivalent widths in observers frame and rest frame, the line flux, the blueshift-adjusted rest wavelengths 
from Tytler & Fan (1992), and the corresponding redshift. The error associated with the placement of the continuum is included 
in the errors on W 0 \, s , TF rest and the line flux. The Lya line is heavily absorbed, so the redshift derived from this line (shown in 
square brackets) is not used in the calculation of the systemic QSO redshift. 


Line 

Arest 

(A) 

A vac 

(A) 

Webs 

(A) 

Wrest 

(A) 

Line flux 

(to -16 erg s 1 cm -2 

Atf 

) (A) 

Ztf 

Lya 

1215.67 

4918 ± 20 

281 ± 15 

69 ±4 

541 ± 33 

1214.97 ± 0.07 

[3.048 ± 0.016] 

Nv 

1240.15 

5005 ± 5 

96 ±4 

24 ± 1.0 

183 ± 27 

1239.16 ±0.28 

3.039 ± 0.004 

Sill 

1263.31 

5101 ±7 

6 ± 2 

1.5 ±0.5 

10 ±3 

— 

— 

Sill/Ol 

1305.57 

5281 ± 7 

15 ±2 

3.8 ±0.6 

22.7 ± 1.5 

1304.24 ± 0.30 

3.049 ± 0.005 

Cii 

1335.31 

5400 ± 3 

9 ± 2 

2.2 ±0.5 

11 ±2 

— 

— 

Si IV /O IV] 

1399.55 

5647 ± 3 

62 ±9 

15 ±2 

81 ±4 

1398.62 ±0.19 

3.037 ±0.002 

Civ 

1549.05 

6249 ± 5 

180 ± 15 

44 ±4 

195 ± 11 

1547.46 ± 0.04 

3.038 ± 0.003 

C hi] 

1908.73 

7722 ± 7 

100 ± 12 

25 ±3 

94 ±9 

1906.53 ±0.09 

3.050 ± 0.004 


a Taken from Wilkes (2000). 



Fig. 8. Zoom on the normalized spectrum of Q1205-30 between the Lya and C IV QSO emission lines. The numbering of the 
absorption lines correspond to the numbering in Tabled 


of galaxies surrounded by a large, homogeneous gaseous en¬ 
velope which is responsible for the absorption, but doubt has 
arisen whether the galaxies at large projected distances are the 
true absorbers (Yanny & York 1992; Vreeswijk et al. 2003; 
Jakobsson et al. 2004). 

We have looked for absorption lines in the QSO spectrum 
due to the foreground galaxies gl and g2 located at impact pa¬ 
rameters 12.6 kpc and 21.3 kpc, respectively. Considering the 
redshifts of gl and g2, the only lines covered by our spectrum 
are Call H and K, Mg II AA2796, 2804, Mg I A2853 for both 
galaxies, as well as Fe II AA2344, 2374, 2383, 2587, 2600 for 
g2. Lines in the Lyman-forest are not well-constrained due to 
the large number of Lya lines. We fitted and removed lines 
corresponding to Lya, Ly (3 and O VI absorption from the C IV 
systems listed in Table 0 In the residuals we identify two 
absorption lines (see Fig. m which may be due to Mg II at 
2 = 0.4722. Alternatively, they could be two Lyman-forest 
lines. In Table [6] we list the 2cr upper limits on the equivalent 
widths of the other lines. 


We find no absorption lines in the QSO spectrum outside 
the Lyman-forest due to the two foreground galaxies gl and g2 
down to a 2cr upper limit on the equivalent widths of ~ 90 mA 
(rest frame), but we cannot exclude that gl is a Mg II absorber 
with equivalent width less than ~ 1 A, as the lines would fall in 
the Lyman forest. Assuming that the g 1 Mg II doublet candidate 
is correctly identified, we cannot exclude that gl is a DLA (see 
Fig. 24 in Rao & Turnshek2000). Conversely, it is unlikely that 
g2 is a DLA system (Fig. 25 and 26 in Rao & Turnshek 2000). 


Old elliptical galaxies have recently been found up to a red¬ 
shift 2 = 1.9 (Cimatti et al. 2004), but not much is known 
about the gas content in elliptical galaxies at these early times. 
The combination of a small impact parameter (corresponding 
to 21.3 kpc) and strict upper limits on absorption for the ellip¬ 
tical galaxy g2 indicates that it has no significant gas column 
density at these radii. The redshift 2 = 0.865 of g2 corresponds 
to a lookback time of 7.1 Gyr. 
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Table 4. Identified absorption lines redwards of Ly a. The 
columns show the line number, the observed wavelength in vac¬ 
uum, the observed equivalent width, the inferred redshift and 
the absorption line and system identification. The error bars 
are ler errors due to read-out noise and photon statistics only. 
The error associated with the position of the continuum is not 
considered. The line numbers refer to the numbering in Fig. [8] 


# 

Avac 

(A) 

w obs 

(A) 

■^abs 

Identification 


1 

4945.0 

1.92 ± 0.02 

2.9917 

Nv 

I 

2 

4963.0 

3.82 ± 0.02 

2.9934/3.0063 

Nv 

I/J 

3 

4979.4 

2.16 ± 0.02 

3.0066 

Nv 

J 

4 

4989.6 

0.18 ± 0.02 

2.8318 

Oi 

E 

5 

4996.8 

0.03 ± 0.02 

2.8308 

Sill 

E 

6 

5067.6 

0.83 ± 0.03 

2.2732 

Civ 

A 

7 

5075.7 

0.57 ± 0.02 

2.2730 

Civ 

A 

8 

5083.5 

1.02 ± 0.03 

3.0332 

Sill 

K 

9 

5110.4 

0.20 ± 0.02 

2.8293 

ClI 

E 

10 

5118.5 

0.73 ± 0.03 

2.6725 

Si IV 

D 

11 

5134.1 

0.20 ± 0.03 

2.9427 

Oi 

G 

12 

5142.2 

0.10 ± 0.02 

2.9423 

Sill 

G 

13 

5152.5 

0.50 ± 0.04 

2.6731 

Si IV 

D 

14 

5187.2 

0.07 ± 0.03 

2.8869 

Cll 

F 

15 

5195.0 

0.04 ± 0.02 

2.9895 

Oi 

I 

16 

5201.9 

0.01 ± 0.02 

2.9880 

Si 11 

I 

17 

5215.3 

0.06 ± 0.03 

3.0051 

Oi 

J 

18 

5251.8 

0.40 ± 0.03 

3.0331 

Oi 

K 

19 

5261.6 

0.05 ± 0.03 

2.9426/3.0338 

CII / Si 11 

G/K 

20 

5338.7 

0.23 ± 0.03 

2.8305 

Si IV 

E 

21 

5347.0 

1.48 ± 0.03 

2.4537 

Civ 

B 

22 

5356.1 

0.97 ± 0.03 

2.4538 

Civ 

B 

23 

5373.3 

0.13 ± 0.03 

2.8305 

Si iv 

E 

24 

5382.5 

0.97 ± 0.03 

3.0332 

Cll 

K 

25 

5414.5 

0.04 ± 0.02 

2.8848 

Si IV 

F 

26 

5448.5 

0.07 ± 0.03 

2.8841 

Si IV 

F 

27 

5495.2 

0.34 ± 0.03 

2.9427 

Si IV 

G 

28 

5531.5 

0.29 ± 0.03 

2.9433 

Si IV 

G 

29 

5563.2 

0.22 ± 0.03 

2.9915 

Si IV 

I 

30 

5577.1 

0.28 ± 0.05 

2.6023 

Civ 

c 

31 

5587.2 

0.13 ± 0.02 

2.6028 

Civ 

c 

32 

5621.7 

1.05 ± 0.03 

3.0335 

Si iv 

K 

33 

5657.9 

0.63 ± 0.02 

3.0334 

Si iv 

K 

34 

5686.9 

2.29 ± 0.02 

2.6732 

Civ 

D 

35 

5696.4 

1.41 ± 0.02 

2.6732 

Civ 

D 

36 

5932.1 

0.58 ± 0.02 

2.8316 

Civ 

E 

37 

5942.5 

0.43 ± 0.03 

2.8320 

Civ 

E 

38 

6015.8 

0.92 ± 0.02 

2.8857 

Civ 

F 

39 

6025.8 

0.49 ± 0.02 

2.8857 

Civ 

F 

40 

6103.9 

1.69 ± 0.02 

2.9425 

Civ 

G 

41 

6114.8 

1.12 ± 0.02 

2.9430 

Civ 

G 

42 

6137.3 

1.19 ± 0.02 

2.9641 

Civ 

H 

43 

6146.9 

1.01 ± 0.33 

2.9637 

Civ 

H 

44 

6159.3 

0.36 ± 0.02 

3.0344 

Si II 

K 

45 

6165.9 

0.10 ± 0.02 

2.8334 

Fell 

E 

46 

6180.0 

1.41 ± 0.02 

2.9918 

Civ 

I 

47 

6190.2 

0.97 ± 0.02 

2.9917 

Civ 

I 

48 

6205.3 

2.52 ± 0.02 

3.0081 

Civ 

J 

49 

6215.7 

2.20 ± 0.02 

3.0081 

Civ 

J 

50 

6245.2 

1.55 ± 0.01 

3.0338 

Civ 

K 

51 

6256.0 

1.47 ± 0.01 

3.0341 

Civ 

K 

52 

6267.3 

0.29 ± 0.01 

3.0481 

Civ 

L 

53 

6278.8 

0.08 ± 0.01 

3.0488 

Civ 

L 


Table 5. Overview of the absorption systems in the spectrum 
of Q1205-30. 


System 

Detected lines 

•Zabs 

A 

Civ AA1548, 1550 

2.2731 ± 0.0001 

B 

Civ AA1548, 1550 

2.4538 ± 0.0001 

C 

Civ AA1548, 1550 

2.6026 ± 0.0004 

D 

Si iv AA1393,1402, Civ AA1548, 
1550 

2.6730 ± 0.0004 

E 

Oi A1302, Sin A1304, Cll 
A1334, Si iv AA1393, 1402, Civ 
AA1548, 1550, Fell A1608 

2.8309 ± 0.0009 

F 

Cll A1334, Si iv AA1393, 1402, 
Civ AA1548, 1550 

2.885 ±0.001 

G 

Oi A1302, Sin A1304, Cll 
A1334, Si iv AA1393, 1402, Civ 
AA1548, 1550 

2.9427 ±0.0003 

H 

Civ AA1548, 1550 

2.9639 ± 0.0003 

I 

Nv AA1238, 1242, Oi A1302, 
Sin A1304, Si iv A1393“, Civ 
AA1548, 1550 

2.991 ±0.002 

J 

Oi A1302, Nv AA1238, 1242, 
Civ AA1548, 1550 

3.007 ±0.001 

K 

Sin AA1260, 1304, 1526, Oi 
A1302, C II A1334, Siiv AA1393, 
1402, Civ AA1548, 1550 

3.0336 ± 0.0004 

L 

Civ AA1548,1550 

3.0485 ± 0.0005 


a Si IV A1403 was too weak to be detected. 


Table 6. Absorption line search in the QSO spectrum due to gl 
and g2. The quoted upper limits are 2<r limits on the equivalent 
widths. The error associated with the placement of the con¬ 
tinuum is taken into account for lines not in the Lyman-forest 
(A obs > 4900 A). 



Ion 

Arest 

(A) 

Aexpected 

(A) 

Wrest 

(A) 

gl 

Mg II 

2796.4 

4116.8 

1.2“ 


Mg II 

2803.5 

4127.4 

1.1“ 


Mg I 

2853.0 

4200.1 

< 0.20 


Call 

3934.7 

5792.8 

< 0.095 


Call 

3969.6 

5844.0 

< 0.068 

g2 

Fell 

2344.2 

4372.0 

< 0.4 


Fell 

2374.4 

4428.4 

< 0.4 


Fe II 

2382.8 

4443.9 

< 0.5 


Fe II 

2586.6 

4824.1 

< 0.5 


Fe II 

2600.2 

4849.3 

< 0.6 


Mg II 

2796.4 

5215.2 

< 0.097 


Mg II 

2803.5 

5228.6 

< 0.107 


Mg I 

2853.0 

5320.8 

< 0.059 


Call 

3934.7 

7338.4 

< 0.089 


Call 

3969.6 

7403.3 

< 0.073 


a The line is in the Lyman-forest, so it may be Ly a. at an intermediate redshift. The 
identification is therefore not secure. 

4.1. Gravitational tensing 

The position of g2 could introduce a lensing effect on S6. 
Knowing the true redshift of g2, we can repeat the calculation 
in Paper I of the radius of its Einstein ring. For easy compar¬ 
ison to Paper I we assume that g2 has a singular isothermal 
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Fig. 9. Velocities of lines originating from absorption system K 
(the LLS) relative to zqi = 3.0331. The model plotted for Lyo 
and Ly/3 has an Hi column density of log Nh , = 19.5. The 
Si II A1304 transition is not included here due to the probable 
blend with C II A1334 at z = 2.9426 (see TableEJ. 

mass distribution with a velocity dispersion a = 300 km s -1 , 
in which case the radius of the Einstein ring is 

0 E =47r(—) ^^ = 1.4 arcsec (1) 

V c / as 

(from the equation following Eq. 4.14 in Peacock 1999), where 
dss is the lens-source angular distance and (Is is the observer- 



Fig. 10. Spectrum of Q1205-30 normalized to the continuum. 
In the spectrum we have subtracted lines corresponding to Lya, 
Ly/3 and O VI absorption from the systems listed in Table [5] 
The dashed line shows the model of a possible Mg II doublet 
from gl. 

source angular distance. Gravitational lensing will introduce 
stretching and distortion perpendicular to the radius vector 
from g2, most notably at a distance of one Einstein radius. 
Emission appearing within one Einstein radius of the centre 
of g2 (roughly corresponding to the dotted circle in the lower 
right part of Fig.[0 is expected to originate from the same small 
region. This stretching effect would result in a flat velocity pro¬ 
file along PA1, which is seen in Fig. Et over a distance of up to 
4 arcsec from the QSO centre towards g2. 

The fact that we only detect one lensed image of the quasar 
constrains the mass of g2 within a projected distance corre¬ 
sponding to the g2 - quasar angular distance (2.77 arcsec). 
Following the calculations of Le Brun et al. (2000), we may 
calculate a model independent upper limit on the projected 
mass M,^ ax (< n) enclosed within a radius n = d E 0\ = 21.3 
kpc: 

ML(< n) = 7TrfEcrft = 2.8 x 1O 12 M 0 , (2) 

where E cr i t = ifc d^d LS * s t ^ le cr ^i ca ^ surface density. This 
corresponds to an upper limit to the radius of the Einstein ring 
of 0 E < 2.77 arcsec. 

Conversely, the following lower limit to the size of the 
Einstein ring allows us to constrain from below the projected 
mass. From Fig. Et it is evident that the Einstein radius is at 
least 0.8 arcsec (the constant part of the velocity profile be¬ 
tween galaxy g2 at ~ 2.8 arcsec and out to the most distant 
measurement at 3.6 arcsec). Thus 0 E > 0 m i n = 0.8 arcsec, 
and 

M{< n) > M(< d L 9 E ) > ^ c2 dsdL = 2.3 x 1O U M 0 . 

4G cLls 

(3) 

Thus the projected mass of galaxy g2 is between 2.3 x 
10 11 Mq and 2.8 x 1O 12 M 0 within a radius of 2.77 arcsec (21.3 

























10 


M. Weidinger et al.: The extended Lyman-a emission around the radio-quiet QS01205-30 


kpc). This is consistent with the super -Mg finding of Sect. 0 
In terms of velocity dispersion the range is 230 km s” 1 < cr < 
420 km s’ 1 , assuming a singular isothermal mass distribution. 


5. A model of the extended Ly a emission 

In trying to understand the extended Ly a emission, we have 
constructed a numerical model, where the quasar lies in the 
centre of a large, optically thin Hi cloud. This model was al¬ 
ready prosposed in Weidinger et al. (2004; hereafter Paper II) 
together with the main conclusions. The details of our calcu¬ 
lations were not included in that paper, but they will now be 
presented here. 

The quasar emission is collimated in a cone with a full 
opening angle T, and the system is seen under an inclination 
angle 9 (see Fig. 3a of Paper II). The HI within the cone is 
photoionized by the quasar UV photons, causing it to emit Lya 
photons when recombining. For the calculation of the extended 
Lya-emission, we follow the treatment of the optically thin 
case given in Gould & Weinberg (1996). The ionization rate 
of a hydrogen atom at a distance r from the quasar is 

pOO 

r(r) = / (t>{v)a(v)dv , (4) 

J VL 

where is the photon flux density 

and the HI ionization cross section is 

/ \ _2 ' 75 

a(u) « 6.3 x HT 18 cm 2 ( — j . (6) 

Here hv l = 13.6 eV is the hydrogen ionization potential, and 
/„(r) is the frequency-specific flux at a distance r from the 
QSO. The production rate of Lya photons per unit volume is 

tlLya(f) = t?thintlHi(r-)r(r), (7) 


where // thin = 0.42 is the fraction of recombinations that result 
in a Lya photon. We assume a power-law hydrogen density 

profile, riHi(r) = n H 1,1 (i^) ■ By integrating along the 

line of sight, /, we obtain the surface brightness (in erg s’ 1 
cm’ 2 arcsec’ 2 ) 


E 


Lya — 


J -^LyatiLya (‘r'jdl D~^ 

47 tD 2 l ~dtt' 


( 8 ) 


Here E^ ya = 10.4 eV is the energy of a Lya photon, and 
Dl the quasar luminosity distance. Da is the angular distance, 
such that jS- is the conversion from cm 2 to arcsec 2 . 

The flux emitted from the quasar close to the Lyman-limit 
frequency, is given by 


fA r ) 



inside the cone, 
outside the cone. 


(9) 


We obtain the observed Lyman-limit flux, f VL)0 bs> using the 
slope of the quasar spectrum and an observed flux, /„ 0i 0 b s , on 
the continuum at z/q, that is 


fv L , ob* = Uo,obs ■ ( 10 ) 

The free parameters of the model are the observed flux at 
i/ 0 , fvo, obs> the spectral slope, a, the redshift, z, the angles 
T and 9, the neutral hydrogen density at a distance of 1 kpc, 
riHi,i, and the density slope, 7 . 


5.1. Model parameters 

Several free parameters are determined directly via observa¬ 
tions. We measure the observed continuum flux at a given point 
in the spectrum. For the spectral index a we use the value found 
in the flux calibration in Sect. 12.31 The obtained values for these 
parameters are listed in Table Q In Paper II a method to obtain 
a relation between the opening angle and the inclination angle 
is described. The free parameters in our model are now reduced 
to the opening angle, T, the hydrogen density scale, «hi.i, and 
the density slope, 7 . 

The numerical implementation of the model was carried out 
by incorporating the cone into a cubic grid inclined at an angle 
9. The Lya production rate was calculated in each grid point 
within the cone, and the process of integrating along the line of 
sight was simply to sum along the grid z-axis. 


Table 7. Input parameters of the model. 


Description 

Symbol 

Value 

Observed flux 

f Aq, obs 

1.3 X HT 16 ergs -1 cm -2 A -1 

Reference wavelength 

Ao 

5946 A 

Spectral slope 

a 

0.65 

Redshift 

z 

3.041 


Table 8. Best-fit values for 9, iihi,l 7, and M 200 for various 
opening angles. 


T 

e 

hhi,i 

( 10 -8 cm -3 ) 

7 

M 200 

(1O 12 M 0 ) 

90° 

22.5° 

9.8 

0.018 

2.6 

110° 

27.2° 

9.3 

0.072 

3.5 

O 

O 

CO 

31.8° 

8.2 

0.093 

5.3 

0 

0 

10 

36.5° 

6.4 

0.064 

6.5 

170° 

41.1° 

5.5 

0.054 

6.5 


5.2. Results of the model 

Using the parameters listed in Table Q and assuming a given 
opening angle and corresponding best-fit inclination angle, we 
now fit the model to the observed surface brightness profile, 
leaving the hydrogen density scale, hhi,i, and slope, 7 , as free 
parameters. The best-fit values for each opening angle are listed 
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arcsec 


Fig. 11. Best-fit surface brightness profiles as seen through a 

I. 2 arcsec wide slit with a seeing of 0.7 arcsec for opening 
angles of 90° (dotted), 110° (solid), 130° (long dashes), 150° 
(solid), and 170° (short dashes). The model parameters are 
listed in Table [8] Overplotted is the measured surface bright¬ 
ness profile of S6 shown in Fig. Eb. Note that three marginal 
measurements at negative offsets are included. 

in Table[8] and the corresponding surface brightness profiles are 
plotted in Fig.m 

The numerical model described in Paper II and above is 
purely geometrical. Assuming in addition that the gas is in free 
fall into a dark matter (DM) halo introduces a velocity field in 
the gas. In Paper II this was used to fit the virial mass, M 200 , 
to the observed velocity profile. The obtained virial masses are 
also listed in Table [8] The conclusions are presented in Paper 

II. 

6. Discussion 

Returning to the main topic of the origin of the extended Lya 
emission, there are several elements of this puzzle that now 
can to be pieced together. There is a LLS at z = 3.0336 ± 
0.0004 (absorber K, see Table an additional absorber at 
z = 3.039 ± 0.002 (Kl, see Sect. 13.3f . the systemic redshift 
of the quasar is z = 3.041 ± 0.001 (Sect. the absorber 
L at z = 3.0485 ± 0.0005 (Table |5j and finally the extended 
Lya emission S6 at z ss 3.049 (see caption of Fig.Et-d). We 
assume that the high-ionization, line locked systems I and J are 
intrinsic to the central engine and will not consider them here. 

6.1. The origin of the extended Lya emission 

The separation of the low- and high-ionization lines in the LLS 
suggests that it is heated by the quasar UV photons and ex¬ 
panding. However, the large velocity offset between the LLS 
and the extended Lya emission (~ 1000 km s -1 ) makes it 
very unlikely that the two are physically associated as previ¬ 
ously thought (Paper I). Since the LLS absorbs all UV photons, 
the extended emission must be located between the LLS and 
the QSO, which makes an association between the extended 


Lya emission and the QSO the most likely one. It was sug¬ 
gested by Haiman & Rees (2001) that neutral gas falling into 
the dark matter halo around a quasar could be photoionized by 
the quasar UV flux, causing the gas to emit Lya photons. 

Interpreting the observations in the Haiman & Rees (2001) 
picture as described in the numerical model (Sect. |5j Paper II), 
we identify Kl as absorption due to a large hydrogen cloud 
around the quasar, in which case the redshift of Kl should be 
close to that of the quasar. A part of this large cloud is pulled 
into the DM halo of the QSO where it is photoionized and gives 
rise to the Lya emission of S6. S6 therefore lies between us and 
the quasar, and the higher redshift is due to infall. The velocity 
of S6 is ~ 500 km s -1 relative to the quasar and the surround¬ 
ing cloud (Kl). The absorption system L is located somewhere 
between the emitting and the absorbing part of the hydrogen 
cloud. There are at least two possibilities for placing the LLS. 
i) It is very close to the quasar and moving at high velocity 
(~ —500 km s -1 ). In this case the LLS has to be very small 
in order for the majority of the UV photons to pass by and 
photoionize the surrounding hydrogen cloud, ii) The LLS is 
sufficiently distant for the Lya photons produced in the cloud 
surrounding the quasar to be redshifted out of the resonance 
wavelength and pass through the LLS unhindered. Because the 
low- and high-ionization lines in the LLS are only mildly sep¬ 
arated in velocity space (~ 60 km s 1 ) we favour the latter 
possibility. 

The extended Lya emission could possibly be explained 
by other scenarios rather than the projected ionization cone. 
However, the combination of imaging and spectroscopic data 
enables us to rule out most of these other scenarios (Paper II). 
i) Jets are believed to be present in radio-quiet quasars. They 
are predicted to extend out to only ~ 0.1 kpc (Blundell et 
al. 2003), more than two orders of magnitude less than the ~ 30 
kpc extent of the emission around Q1205-30. ii) Outflowing 
galactic winds are generally thought to be triggered by the 
cumulative effect of many supernovae exploding inside the 
galaxy, which would metal-enrich the outflowing gas. Around 
radio-loud quasars this is typically seen as extended CIV emis¬ 
sion with strength 7 — 10% of that of the extended Lya line 
(Heckman et al. 1991b). Our detection limit at the position 
of redshifted Civ is 4 x 10 -18 erg s -1 cm -2 arcsec -2 (3cr), 
which would have enabled us to detect the typical Civ line 
seen around some radio-loud quasars (a strength of 10% of the 
Lya line corresponds to 7 x 10 -18 erg s -1 cm -2 arcsec -2 ). 
The fact that we do not detect this line makes it unlikely that 
it is a supernova powered outflow. The most plausible explana¬ 
tion that remains is cosmological infall of hydrogen. 

6.2. Extending the sample 

To date only a handful of detections of extended Lya emission 
around RQQs has been reported (Steidel et al. 1991; Bremer 
et al. 1992; Hu et al. 1996; Petitjean et al. 1996; Bunker et 
al. 2003). Hu et al. (1991) find no extended Lya emission in 
their sample of seven radio-quiet quasars down to a limiting 
flux of 2 x 10 -16 erg s -1 cm -2 . We have compiled a list 
of fluxes and luminosities for the detections of extended Lya 
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Table 9. Redshifts, fluxes and luminosities of the Lya emission from the host galaxies of radio-quiet quasars. 


Quasar 

z 

Lya flux 

(10 -17 erg s 1 cm -2 ) 

Lya luminosity 

(10 42 ergs -1 ) 

Reference 

Q 0054-284 

3.616 

6 

7.3 

Bremer et al. (1992) 

Q 0055-26 

3.656 

8 

10 

Bremer et al. (1992) 

Q1548+0917 

2.746 

47 

29 

Steidel et al. (1991) 

Q 1205-30 

3.041 

70 

56 

This work 

BR 1202-0725 

4.69 

26 

58 

Hu et al. (1996), Petitjean et al. (1996) 


emission (see Table [h}. We measure an average Lya surface 
brightness of ~ 7 x 10 -1 ' erg s -1 cm -2 arcsec -2 around 
Q1205-30. Assuming a spatial extent of 10 arcsec 2 , we infer 
a Lya flux of ~ 7 x 10 -16 erg s -1 cm -2 . The luminosity in 
the Lya line of S6 is at the high end of the ones listed in Table 

El 

It is important for the study of the link between galaxy 
and quasar formation to understand how frequent and under 
which circumstances extended Lya emission arises around 
RQQs. The most efficient method to detect Lya emission 
around QSOs is narrow band imaging. The study of ex¬ 
tended Lya emission can successfully be combined with nar¬ 
row band searches for Lya emitting (proto)-galaxies around 
QSOs (e.g. Fynbo et al. 2001; Fynbo et al. 2003). The follow¬ 
up multi-object spectroscopy needed to confirm candidate Lya 
emitters may conveniently be utilized to study any extended 
Lya emission associated with the QSO. Alternatively, integral 
field spectroscopy is a promising method for very detailed stud¬ 
ies of extended emission at high redshifts (Bower et al. 2004). 
The method makes it possible to map out the entire veloc¬ 
ity field of the extended emission, strongly constraining any 
model. 

It is imperative for any search for extended Lya emission 
around RQQs to go to very deep detection limits. Haiman & 
Rees (2001) predict that haloes of infalling gas around quasars 
should be seen in Lya emission with a typical surface bright¬ 
ness around 10 -18 — 10 —1 ' erg s -1 cm -2 arcsec -2 and typi¬ 
cal angular sizes between 2 — 3 arcsec, i.e. with typical fluxes 
around 3 x 10 -18 — 7 x 10 -17 erg s -1 cm -2 . This limit has 
only been reached for very few surveys. 

A larger sample of QSOs with extended Lya emission will 
make it possible to address primary issues like morphology, 
environment, luminosity function etc. The corresponding DM 
halo masses obtained in a similar fashion as in Paper II may 
be compared to black hole masses obtained via the correlation 
between Mbh and quasar emission line widths (Vestergaard 
2002; McLure & Jarvis 2002). A Mdm — Mbh correlation 
could provide a powerful consistency check of iV-body hydro- 
dynamical simulations. 

We are currently looking for extended Lya emission in a 
small sample of quasars using an analysis similar to what has 
been employed in this paper. 
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